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Проведен анализ следящего оптического фазометра с временным дискриминатором видеосигнала, для которого получены оп- 
тимальные соотношения параметров при обработке видеосигнала малоконтрастного интерференционного поля, то есть, в усло- 
виях гауссова распределения и марковской модели сигнала и фона. Для пуассоновской модели распределения сигнала и фона 
приведены выражения, определяющие оптимальные размер и скорость следящего сканирования малоконтрастной интерфе- 
ренционной картины по условию минимума погрешности слежения за полосой интерференции. 


Средства интерферометрии находят широкое при- 
менение в ряде оптических и радиотехнических уст- 
ройств и систем [1-7]. Интерес к интерферометричес- 
ким средствам, обеспечивающим прием малоконтра- 
стных низкой интенсивности изображений в услови- 
ях воздействия помех, особенно возрос в последние 
годы в связи с появлением персональных компьюте- 
ров высокой мощности и новых методов, например 
[2, 3, 5, 6], вейвлет-анализа сигналов и изображений, 
делающих возможным обработку картин интерфе- 
ренции со скоростью, достаточной для получения в 
реальном времени информации о динамике измене- 
ния трехмерной структуры фазового распределения 
изображений наблюдаемого средством пространства 
[2^1]. Задача получения оценок динамики изменения 
пространственно-временных характеристик волново- 
го фронта оптического излучения значительно слож- 
нее подобной [5, 6] радиоволновой задачи. Это объяс- 
няется тем, что фотоприемные (ФП) устройства опти- 
ческого диапазона не осуществляют непосредственно 
регистрацию комплексной амплитуды волнового по- 
ля. Известные методы регистрации комплексной 
амплитуды, определяя способы получения простран- 
ственных и временных характеристик оптического 
поля, являются косвенными и основаны [1, 2, 4] на 
обширной гамме различных модификаций основного 
принципа интерференционного взаимодействия от- 
дельных частей когерентного или некогерентного оп- 
тического поля. Название методов, например, голог- 
рафическая, двухэкспозиционная, спекл-, сдвиговая, 
муаровая, гетеродинная, интерферометрия интенсив- 
ностей отражают особенности аналитических подхо- 
дов, дающих возможность регистрации параметров 
исследуемого явления, с точностью и разрешением, 
определенными размерностью длины волны исполь- 
зуемого излучения. Так, интерферометрические уст- 
ройства контроля формы поверхностей по измерению 
фазы фронта потенциально позволяют определять 
комплексную амплитуду оптического поля и другие 
информативные параметры как апериодических, так 
и колебательных сверхмалых ~(10 2 ...10 3 )Я перемеще- 
ний поверхности [4]. В этом случае изменения фазо- 
вого фронта наблюдаются как результат сложения 
волны исследуемого поля с полем другого опорного 
когерентного оптического излучения. Регистрация 
этих изменений фазы фронта заключается в последу- 
ющем изучении структуры, контраста, динамики по- 


ведения и положения интерференционных полос 
(ИП) в плоскости изображения. Малый контраст 
изображения ИП при слабых световых сигналах явля- 
ется одним из основных факторов ограничения эф- 
фективности средств, учитываемый при проектиро- 
вании систем адаптивной оптики и фазометрических 
систем контроля в реальном времени динамической 
деформации волнового фронта. 

Эта проблема порождает необходимость удов- 
летворения противоречивых требований к фазо- 
метрической системе, а именно, обеспечения мак- 
симальной эффективности фильтрации помех при 
высоких динамической и статической точности в 
измерении и воспроизведения фазы волнового 
фронта анализируемого излучения. Наиболее эф- 
фективное решение этой проблемы возможно 
только на основе синтеза системы, реализующей 
минимум среднего квадрата случайной погрешнос- 
ти, обусловленной воздействием искажающих 
факторов - помех, при одновременном удовлетво- 
рении технических требований к качеству ее других 
характеристик, например, обеспечения наиболь- 
шей чувствительности. Вместе с тем, вопросы син- 
теза оптических фазометров с определенными ха- 
рактеристиками практически не разработаны и в 
литературе недостаточно освещены. Основные ре- 
зультаты, полученные в этом направлении и изло- 
женные, например, в [1], не отражают всей полно- 
ты требований, предъявляемых на практике к по- 
казателям качества таких систем. 

Задача синтеза оптимального оптического фа- 
зометра для контроля параметров волнового фрон- 
та излучения с математической точки зрения мо- 
жет быть рассмотрена как экстремальная задача. 
Она сводится к исследованию на условный мини- 
мум функционала, определяющего средний квад- 
рат случайной погрешности (расстройки или не- 
вязки) системы при некоторых ограничивающих 
условиях, которые соответствуют требованиям 
чувствительности, статической и динамической 
точности оптического фазометра. Так, для рас- 
смотренных в [8, 9] следящих фазометров размер и 
скорость следящего сканирования могут быть оп- 
тимизированы на основе общепринятой теории 
статистических решений по критерию минимума 
средней квадратической погрешности, максимизи- 
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рующей точность системы слежения и зависящей 
от энергетических, например, отношения сиг- 
нал/шум и неэнергетических ее параметров. Рас- 
смотрим их влияние подробнее путем исследова- 
ния эффективности применения в фазометрах 
принципа временного дискриминирования совме- 
стно с интерференционным методом контроля фа- 
зы волнового поля. В фазометр этого типа заложен 
принцип слежения за подвижным объектом [10] с 
время-импульсным способом формирования сиг- 
нала ошибки [11], основанного на том, что смеще- 
ние центров следящего сканирования и изображе- 
ния ИП приводит к временному сдвигу видеоим- 
пульса относительно опорного уровня - виртуаль- 
ного центра сканирования. Измеряя временной 
сдвиг видеоимпульса от ИП относительно сигнала 
сканирования для формирования определенного 
величиной этого сдвига аналогового или цифрово- 
го сигнала ошибки положения ИП, можно совмес- 
тить центры сканирования и изображения ИП. 

В представленном ниже варианте слежения за 
смещением ИП, в течение периода или интервала 
следящего сканирования Т ск , ИП пересекает эле- 
мент ФП один раз, и выходной сигнал ФП устрой- 
ства будет представлять собой результат свертки 
распределений плотности потока фотонов в изобра- 
жении ИП размером 2 г 0 и функции чувствительнос- 
ти апертуры ФП, с учетом формы его элемента раз- 
мером 2 г а . Для общности рассматриваем случай 
замкнутой области связной структуры фрактала 
изображения ИП некоторого порядка интерферен- 
ции. При распределении энергетической освещен- 
ности в изображении пятна ИП, в том числе нулево- 
го порядка кольцевой структуры, в плоскости ФП 
Е(х',У) и функции совместного распределения 
прозрачности оптической апертуры и чувствитель- 
ности элемента ФП р а (х,у) амплитуда сигнала на его 
выходе определится двумерным интегралом свертки 

іс = е Е 0 л го а 2 (1+ а 2 у 1 ехр [-(ті + пі) / (I + а 1 )\, (1) 


где г - чувствительность ФП (в общем случае функ- 
ция Я); а=г а /г й , піігХъ/гц, п 0 =у п /г () - относительные 
размер элемента ФП и координаты; Е (І - энергети- 
ческая освещенность сигнала (в точке х 0 =О; у 0 =0), с 
максимальным значением сигнала на выходе ФП 

і с =еЕ 0 лгоа 2 /(1 + а 2 ). (2) 


При равномерной прозрачности апертуры ФП 
круглой формы радиусом г а 


Ра(х,у) = 


\,для\х\< л/ г 2 -у 2 , | у | < уіг 2 -х 2 ; 
О, для | х [ > ^г 2 -у 2 , \у\> V г 2 -х 2 , 


величина сигнала на основе (1) и (2) определяется 
выражением 


Іс.а — & Е 0 П /о 


Ѵя , 


и-по 

х | [ехр(-йоі)] 


Ф(^а 2 - (лоі + щ ) 2 — то ) + 
+Ф(^]а 2 — (и„і + По ) 2 + пь) 


СІП а 


(здесь Пт = (у-у 0 )/г 0 ;Ф(и) = -^\ехр(-1 2 )сІ(- ин- 

Ѵ7Г І 

теграл вероятности) с максимальным значением 
сигнала 


і С м = еЕ 0 кг 0 2 [\-ехр(-а 2 )]. (3) 

Для квадратной формы апертуры ФП со сторо- 
ной 2 г кв равномерной прозрачности 

. . \,при\х\<г кв ,\ у\<г т \ 

Ра(х, У ) = 

О, при | х I > г кв , | у | > г кв , 

сигнал на выходе ФП при а=г кв /г 0 равен 

4 = 0,25еЕ 0 я г 0 2 [Ф(йг + »г 0 ) + Ф(а-/я 0 )] х 
х[Ф(й + Я„) + Ф(йГ-й 0 )], 

с максимальным значением і см = еЕ 0 кго [Ф(а)] 2 . (4) 

Если учтем, что с точностью не хуже 0,7 % 
Ф(+м)*+{1-ехр[(-4/я)и 2 ]} 1/2 , то получим макси- 
мальный сигнал, соответствующий наибольшей 
интенсивности 

ісм ~ еЕ 0 пг 0 2 {\ -ехр[(-4/7г)я 2 ]}. (5) 

Часто представляет интерес среднее значение 
сигнала, учитывающее случайность положения 
изображения ИП между тактами сканирования. Для 
одномерной процедуры сканирования с шагом 2 5 Х 
среднее значение сигнала на выходе ФП составит 
8 Х 

іср = | (6) 

-< 5 * 

где в'(х 0 )=1/(2(5 І ) - равномерная плотность распре- 
деления. Для двумерного точечного с шагом 
25 х у.25 у следящего сканирования среднее значение 
сигнала ФП найдется из 

бу д х 

іср=\ | І с (х 0 ,у 0 )\ѵ{х 0 )ц>{у а )с1х 0 с1у 0 , (7) 

-8у -8 Х 

где Цх 0 ,у 0 )=1/(4(5„(5 Ѵ ) Для линейного изображения 
ИП: >ѵ(х 0 )=1/(2(5Д, \ѵ(у 0 )= 1; 5=г а . Для круглой и 
квадратной форм апертур ФП с гауссовской проз- 
рачностью выражения (6) и (7) соответственно 
принимают вид: 

іср =4*Ф(М(Ѵя72&і)ехр(- ті / л/і + й 2 ), 

Іср. 1 — (^ХІ /(2/^x1 )э 

где Ъ й =5У(гуК+ я 2 ). При фоновой плотности потока 
Е ф , величина которого косвенно учитывает темновую 
составляющую шума, сигнал ФП с апертурой круглой 
формы равномерной или гауссовской прозрачности 
равен: і ф =еЕ ф лгі, а на выходе ФП с апертурой квадрат- 
ной формы равномерной прозрачности: і ф =еЕ ф АЕ кв . 

На практике заданными являются яркость фона Б ф 
и сила света излучения или / 0 - интенсивность интер- 
ференционной картины. При этом освещенность 
изображения полосы в плоскости ФП равна: 
Д=г 0 / 0 /(Д/2г 0 ) 2 а фона Е ф = ятДД Д/2/) 2 , где т„ -коэф- 
фициент пропускания среды и объектива; Д -диаметр 
входного зрачка и/ - фокусное расстояние объектива. 
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При следящем сканировании х„(?) ИП (рис. 1 и 2) 
относительно опорного среднего положения 
х м ( Г) =х мо { /,.+ 0 , 5 Т ск ) центра траектории сканирова- 
ния, производимого пьезоприводом типа [12], вы- 
ходной сигнал устройства представляет собой перио- 
дическую последовательность видеоимпульсов. Их 
огибающая промодулирована случайным образом по 
ширине и временному положению видеоимпульсов 
с внутренним заполнением случайным марковским 
потоком однофотонных импульсов, распределение 
интенсивности которых отражает изменение плот- 
ности потока энергии в сканируемом поле ИП. Для 


упрощения формализации анализа точностных ха- 
рактеристик, примем, что для большинства ФП, ис- 
пользуемых для регистрации интерференции, спра- 
ведлива пуассоновская модель сигнала. Тогда, вы- 
ходной сигнал ФП представляет собой пуассоновс- 
кую последовательность наносекундных импульсов 
одноэлектронной природы с плотностью распреде- 
ления вероятности интервала их следования А 1 
рШ) = (1/ г„„)ехр(-Д// Топ), 

где т„„=1/й„(?) - средний период следования выход- 
ных импульсов, величина которого обратно пропор- 



Рис. 1. Структурная схема формирования сигнала рассогласования ВД 




Рис. 2. Динамика переходных процессов при работе следящего фазометра с ВД 
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циональна их интенсивности п и (і), определяемой 
плотностью светового потока в изображении поло- 
сы на чувствительной поверхности ФП. При слеже- 
нии за изображением ИП достаточно большой осве- 
щенности, сигнал ФП имеет вид последовательнос- 
ти видеоимпульсов. Форма их огибающей определя- 
ется видом распределений плотности светового по- 
тока в изображении ИП, прозрачности оптической 
апертуры и элемента ФП и их геометрией. Для соот- 
ветствия принятой для анализа пуассоновской мо- 
дели, будем считать, что полоса пропускания опти- 
мального фильтра (рис. 1) блока обработки выход- 
ного сигнала ФП, в необходимой для реализации 
модели мере, согласована с частотой следования вы- 
ходных импульсов. Среднее значение интенсивнос- 
ти следования импульсов видеосигнала одноэлект- 
ронной природы, как математическое ожидание пу- 
ассоновской случайно распределенной величины, 
равно сумме средних интенсивностей темнового 
шума внешнего фона п ф (і) и сигнального 
п с {і) потоков импульсов: п(1)=Пт(і)+Пф(і)+п с (і), опре- 
деляемых принятыми (аддитивной, мультиплика- 
тивной, аппликативной) моделями изображения. 
Здесь йф(/)=ь-С ф 5ф„/(у и п с (1)=еЕ 0 5 , ф ,/^ вычисляются с 
учетом особенностей сканирования, ранее рассмот- 
ренных [8, 10]; = 1 ,6- 10 19 Кл - заряд электрона; 

5ф„ - площадь элемента ФП. 

Сигнал временного дискриминатора (ВД) на 
интервале т д вычисления величины рассогласова- 
ния положения видеосигнала и центра следящего 
сканирования, отражающий расстройку системы 
или ее разбаланс, представляется в виде [10] 

х д 

ѵ а = к д \{РМ)-Р м (і+тШУи 

О 

где К д - коэффициент пропорциональности, влияю- 
щий на крутизну характеристики преобразования 
ВД; Р т ((,т) - нормированная знакопеременная 
функция, описывающая алгоритм работы ключевых 
каскадов ВД. Она принимает значение "1" при замк- 
нутом состоянии ключа и "0" - при разомкнутом сос- 
тоянии. В рассматриваемом случае, учитывая диск- 
ретные свойства пуассоновского потока импульсов, 
с достаточной степенью точности анализа операцию 
интегрирования можно заменить суммированием на 
том же временном интервале усреднения. Тогда 
среднее значение сигнала рассогласования положе- 
ния видеосигнала относительно центра следящего 
сканирования, определяемое как математическое 
ожидание, и его дисперсия найдутся по формулам 

от„ =К а ^ур і {1)-Р м {і + т)\п(і)М, 

а 2 и =КІ^[Щі) + Р І+1 (< +*№«)&<■ 

Здесь учтены свойства знакопеременной функ- 
ции ІД0;1): /Д/)'/у,(г)=0 и /уП)=Т(/). Когда ско- 
рость следящего сканирования относительно вы- 
сока по сравнению с динамикой смещения отсле- 
живаемой ИП, то величина А г интервала рассогла- 
сования положения видеосигнала и центра траек- 


тории следящего сканирования длиной і ск лежит в 
области малых значений, и выражение для сигнала 
дискриминатора, с использованием методики [11], 
принимает простой вид: т=К д п 0 -т д -Ат/Э ск -Т с1 . Отку- 
да крутизна дискриминатора 

Ѳ„=дІІ/дАт =К д щ-т д /9 ск - Т а 
для т д /Т с =1 составит: Ѳ д =К д щ/ & ск . 

Дисперсия сигнала рассогласования для т д = Т а и 
длительности видеосигнала слежения, определяе- 
мого на уровне е \ т вс =2г^а{\+а)Т сл /І т найдется из 
выражения 

а] = КІ ■Т а ,[Пф'Т +щ(2т{\+а)/ (. ск )\. ( 8 ) 

Нормируя на крутизну ВД среднюю квадрати- 
ческую погрешность рассогласования при следя- 
щем сканировании (время усреднения равно ин- 
тервалу следящего сканирования) а с =ст и /Ѳ д , най- 
дем оптимальное значение скорости и размера сле- 
дящего сканирования из условия минимума пог- 
решности слежения да сл /дЭ с =0: 

&скорІ — Йі й(1 й )й() I Т с , /7ф . 

Здесь й 0 и и® определяются с учетом выражений 
(2-5), отражающих влияние на интенсивность сиг- 
нала формы и вида распределений чувствительнос- 
ти элемента ФП и оптической апертуры. Учитывая 
это для соответствующих вариаций форм апертур 
ФП, получим оптимальный размер траектории 
следящего сканирования І скт = Т ся : 

Р скорі ~ ЙОшах Д) (1 + й) X 

х{[й 2 /(1 + й 2 )];[1 -ехр(-й )];[1 -ехр(Яй/я)]}/йПф тах . 

Поведение полученной из (8) и нормированной 
на значение при а = 1 средней квадратической пог- 
решности для однократного следящего сканирова- 
ния в зависимости от отношения сигнал/фон для 
гауссовых апертуры и элемента ФП и изменении 
соотношения 2 г 0 /і ск орІ , принимающего значения, 
кратные длине волны Я— (4 г 0 /^Я) или в угловой ме- 
ре - кратные л(2ги/І с =2гц/хл), представлено на 
рис. 3. Абсолютное значение средней квадратичес- 
кой временной ошибки слежения за ИП для треу- 
гольной формы сигнала с каниро вания составит: 
Аг а ш =0,25а( 1 +а)(п с /п ф ) 1 Т с ^Т у /2Т а , где Т у - интер- 
вал времени усреднения сигнала ВД, равный [8] дли- 
тельности формирования одной "ступеньки" управ- 
ляющего напряжения основного канала, а отноше- 
ние Т у /2 Т а равно числу усредняемых значений отсче- 
тов сигнала ВД на такте управления пьезоприводом. 

Численное значение верхней границы времен- 
ной погрешности (наихудшее значение при оценке 
точности) метода определим для равных значений 
интенсивности сигнальных и фоновых импульсов 
(й с =й ф ) и элемента пьезопривода [12], выполнен- 
ного на основе пьезопакета ПП-4 0 30x18x5,5 мм 3 
с суммарной емкостью 0,5 мкФ +5 %, резонансной 
частотой / О =40 кГц, допустимым напряжением 
+300 В и чувствительностью 5-10~ 3 мкм/В. Для 
Т у = 10~ 2 с, Т а = 2,5- 10 5 с и а = 1 значение верхней 
границы временной погрешности равно: 


83 




Известия Томского политехнического университета. 2004. Т. 307. № 4 


шрДт„ 0Ш =тг л/2-10 5 с. 

В долях длительности такта интервала управле- 
ния (~10~ 2 с) она составит л: Ѵ2- 1 0 3 или порядка 
0,4 %. С увеличением отношения сигнал/фон по- 
вышается и точность слежения. Для отношения 
сигнал/фон ~10 вносимая шумовыми флуктуация- 
ми погрешность оценки снижается на порядок, и 
верхняя граница временной погрешности слеже- 
ния не превышает 0,04 % длительности одного так- 
та интервала управления пьезоприводом [12]. 



Максимальное отношение сигнат/фон 

Рис. 3. Влияние отношения сигнал/фон и размера скани- 
рования на погрешность оценки координат ПИ 

Основной задачей интерферометрии является 
наиболее полное определение параметров поля ин- 
терференции, формируемого в процессе преобразо- 
вания волнового фронта интерферометром измери- 
тельного средства. Интерференционным картинам 
малого контраста свойственна сложная структура с 
локализованными фрактальными областями, отра- 
жающими свойства объектов и фона в поле наблю- 
дения. Динамика изменения поля наблюдения обус- 
ловливает пространственные и временные вариа- 
ции локализованных фрактальных областей, опре- 
деляющих информационный формализм видоизме- 
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нения общей структуры полного изображения. Это 
определяет необходимость разработки новых мето- 
дов для большей полноты анализа поля интерфе- 
ренции, чтобы определять совокупно с простран- 
ственно-частотными характеристиками поля в це- 
лом [2, 4] его локальные особенности [3] с возмож- 
ностью оценки динамики локальных и глобальных 
вариаций изображения поля интерференции. 

Наиболее широко для анализа поля интерфе- 
ренции используется основанное на габоровском 
представлении аналитических свойств оптических 
сигналов Фурье-преобразование [13]. Для этого 
сигнал задается в области переменных часто неог- 
раниченной протяженности в виде ортогонального 
разложения по частотным базисным функциям. 
Частотные компоненты разложения поля изобра- 
жения определяются путем интегрального преобра- 
зования полной реализации сигнала на бесконеч- 
ном интервале. Значительное усреднение затрудня- 
ет выделение локальных особенностей общей 
структуры поля интерференции. При обработке ма- 
лого контраста картины интерференции зашумлен- 
ные ее локальные вариации практически теряются. 

В настоящее время для анализа, обработки и 
оценки параметров волновых полей идет интенсив- 
ная апробация нового метода [2, 3, 5, 6, 14] многома- 
сштабного анализа сигналов - "вейвлет-анализа". 
Описанный в [8] принцип следящего фазометриро- 
вания с оптимальными размером и скоростью сле- 
дящего сканирования обеспечивает получение как 
локальных фазовых изменений, так и их глобальное 
распределение по полю интерференции при суще- 
ственно сниженном объеме данных, необходимых 
для получения достаточной информации о поле ин- 
терференции. Это позволит осуществлять вейвлет- 
преобразование изображения фазового поля в ре- 
альном времени при меньших вычислительных ре- 
сурсах средств обработки. 
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